Tackling Systematic Errors in Quantum Logic Gates with Composite Rotations 
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We describe the use of composite rotations to combat systematic errors in single qubit quan- 
tum logic gates and discuss three families of composite rotations which can be used to correct 
off-resonance and pulse length errors. Although developed and described within the context of 
NMR quantum computing these sequences should be applicable to any implementation of quantum 
computation. 
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PACS numbers: 03.67.-a, 76.60.-k, 82.56.Jn 



I. INTRODUCTION 

Quantum computers Q| are information processing de- 
vices that use quantum mechanical effects to implement 
algorithms which are not accessible to classical comput- 
ers, and thus to tackle otherwise intractable problems 0. 
Quantum computers are extremely vulnerable to the ef- 
fects of errors, and considerable effort has been expended 
on alleviating the effects of random errors arising from 
decoherence processes (|[ ||[ It is, however, also im- 
portant to consider the effects of systematic errors, which 
arise from reproducible imperfections in the apparatus 
used to implement quantum computations. 

The effects of systematic errors are clearly visible in nu- 
clear magnetic resonance (NMR) experiments |J which 
have been used to implement small quantum computers 
|fl ||, ||, 10, [III 0- Implementing complex quantum al- 
gorithms require a network of many quantum logic gates, 
which for an NMR implementation translates into even 
longer cascades of pulses. In these cases small systematic 
errors in the pulses (which can be ignored in many con- 
ventional NMR experiments) accumulate and have sig- 
nificant effects. 

It makes sense to consider systematic errors as some of 
them can be tackled relatively easily. In the Bloch pic- 
ture, where unitary operations are visualized as rotations 
of the Bloch vector on the unit sphere, systematic errors 
are expressed as rotational imperfections. The sensitiv- 
ity of the final state to these imperfections can be much 
reduced by replacing single rotations with composed ro- 
tations as discussed below. 



II. SYSTEMATIC ERRORS IN NMR 
QUANTUM COMPUTERS 

Any implementation of a quantum computer requires 
quantum bits (qubits) on which the quantum informa- 
tion is stored, and quantum logic gates which act on the 
qubits to process the quantum information. Fortunately 



'Electronic address: jonathan.jones@qubit.org 



it is only necessary to implement a small set of quantum 
logic gates, as more complex operations can be achieved 
by joining these gates together to form logic circuits. A 
simple and convenient set comprises a range of single 
qubit gates together with one or more two qubit gates, 
which implement conditional evolutions and thus logical 
operations (f|. 

NMR quantum computers are implemented |TT| us- 
ing the two spin states of spin-1/2 atomic nuclei in a 
magnetic field as the qubits. Transitions between these 
states, and thus single qubit gates, are achieved by the 
application of radio frequency (RF) pulses. Two qubit 
gates require some sort of spin-spin interaction, which 
in NMR is provided by the scalar spin-spin coupling (J 
coupling) interaction. While this does not have quite the 
form needed for standard two qubit gates, it can be easily 
sculpted into the desired form by combining free evolu- 
tion under the background Hamiltonian (which includes 
spin-spin coupling terms) with the application of single 
qubit gates [ fTD . 

As single qubit gates involve the application of external 
fields they are vulnerable to systematic errors in these 
fields. In the ideal case, the application of a RF field in 
resonance with the corresponding transition with relative 
phase <j> (in the rotating frame ||) will drive the Bloch 
vector through some angle about an axis orthogonal to 
the z-axis and at an angle <f> to the x-axis. The rotation 
angle, 9, depends on the nutation rate induced by the RF 
field, usually written v\, and the duration of the pulse, 
t. In practice the RF field is not ideal, and this leads 
to two important types of systematic errors, pulse length 
errors and off-resonance effects 

Pulse length errors occur when the duration of the RF 
pulse is set incorrectly, or (equivalently) when the RF 
field strength deviates from its nominal value, so that 
the rotation angle achieved deviates from its theoretical 
value. Within NMR this effect is most commonly ob- 
served as a result of spatial inhomogeneity in the applied 
RF field, so that it is impossible for all the spins within a 
macroscopic sample to experience the same rotation an- 
gle. Off-resonance effects arise when the RF field is not 
quite in resonance with the relevant transition, so that 
the rotation occurs around some tilted axis. 

Composite pulses || |lj, are widely used in NMR 
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to minimize the sensitivity of the system to these errors 
by replacing simple rotations with composite rotations 
which are less susceptible to such effects. However, con- 
ventional composite pulse sequences are rarely appropri- 
ate for quantum computation because they usually incor- 
porate assumptions about the initial state of the spins. 
Such starting states are not known for pulses in the mid- 
dle of complex quantum computations, and it is therefore 
necessary to use fully-compensating (type A) composite 
pulse sequences [|l5], which work for any initial state. 
Composite pulses of this kind, which do not offer quite 
the same degree of compensation as is found with more 
conventional sequences, are of little use in conventional 
NMR, and have received relatively little study. They are, 
however, ideally suited to quantum computation. 



III. OFF-RESONANCE ERRORS 

The problem of tackling off-resonance errors was ini- 
tially studied by Tycko [ Em ; his results were then ex- 
tended by Cummins and Jones [|l7| [1^] . Here we describe 
two families of composite pulses which can be used to 
compensate for off-resonance errors, and show how they 
can be derived using quaternions. 

The original method used to develop many type A com- 
posite pulse sequences [I?], [l8| was based on dividing 
the propagator describing the evolution of the quantum 
system into intended and error components, and then 
seeking to minimise the error term. While this approach 
is effective, it is cumbersome, and a much simpler ap- 
proach can be adopted for single qubit gates, which are 
simply rotations on the Bloch sphere and so can be mod- 
elled by quaternions. The quaternions corresponding to 
individual pulses can be multiplied together to give a 
quaternion description of the composite pulse, which can 
then be compared with the quaternion of the ideal sys- 
tem. 

A quaternion is often thought of as a vector with four 
coefficients, but when describing a rotation it is more 
useful to regroup these coefficients as a scalar and a three- 
vector, 

q = {*,v} (f) 

where 

s = cos(0/2) (2) 

depends solely on the rotation angle, 9, and 

v = sin(0/2)a (3) 

depends on both the rotation angle, and a unit vector 
along the rotation axis, a. Thus the quaternion describ- 
ing an on-resonance pulse with phase angle <j) is 

q 9l p = {cos(0/2), sin(0/2){cos(0), sin(0), 0}}. (4) 



An off-resonance pulse is conveniently parameterised by 
its off- resonance fraction / = 8/vx (where 5 is the off- 
resonance frequency, and V\ the nutation rate), and is 
described by the quaternion 

q e0 = {cos(0'/2), !^£2|{cos(0), sin(0), /}} (5) 

where 9' = 9^1 + f 2 , and 9 is now the nominal rotation 
angle, that is the rotation achieved when / = 0. The 
quaternion describing a sequence of pulses is obtained by 
multiplying the quaternions for each pulse according to 
the rule 

q x * q 2 = {si • s 2 - Vi • v 2 , SiV 2 + s 2 v x + Vi A v 2 }. 

(6) 

Finally, two quaternions can be compared using the 
quaternion fidelity |f5| 

-^(qi,q2) = |qi • q2| = |si • s 2 + vi ■ v 2 | (7) 

(it is necessary to take the absolute value, as the two 
quaternions {s,v} and {— s, — v} correspond to equiva- 
lent rotations, differing in their rotation angle by integer 
multiples of 2ir) . 

Following our previous work [jl7| we seek to tackle off- 
resonance errors in a 9 X pulse using a sequence of three 
pulses applied along the x, —x and x axes; pulses with 
any other phase angle can then be trivially derived by 
simply adding the desired value to the phase angles of 
all the pulses in the sequence. Such sequences can be 
described completely by the nominal rotation angles of 
the three pulses, 9\, 9 2 and 9%. The composite quaternion 
for this composite pulse is complicated, but the situation 
can be greatly simplified by expanding it as a Maclaurin 
series in / and neglecting all terms above the first power. 
This gives 




while the ideal quaternion has the form 

{cos(0/2),{sin(0/2),O,O}}. (10) 

ft now remains to chose the three nominal rotation angles 
so that these equations agree. 

First we note that in order to achieve the correct rota- 
tion angle, s — cos(6*/2), we must choose our angles such 
that 9\ — 6*2 + 6*3 = 9 + 2an (where a is any integer) . We 
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also note that the y component of v should equal zero, 
and that this can be achieved by choosing 0j = 8 3 + 2bn 
(where b is any integer). These two choices give 



(11) 



v = <j sin -, 0, / I sin - - 2 sin i - - n L 



Finally we choose 0i such that the z component of v 
equals zero; this gives 



li = arcsm 

1 2 



sin(0/2) \ 
J 



(12) 



Combining this value with our previous relations between 
the angles gives 



. ( zm(6/2) \ 

1\ = ln\-K H arcsm 

2 I 2 / 



6-2 = 2n 2 7r — 2 arcsin 



= 2ri37r + — — arcsin 



(g/2) \ 

sin(0/2)\ 
2 J 



(13) 
(14) 
(15) 



where n 1; n 2 and n% are integers, subject to the phys- 
ical restriction that the resulting pulse angles must be 
positive. 

These solutions have the same general form as those 
found previously fL7|j . Although they appear to differ 
in detail the expressions are, in fact, identical: taking 
the values ni = 1, = 1 and 7J3 — gives our previ- 
ous family of solutions JT^ ], referred to by the acronym 
CORPSE (Compensation for Off-Resonance with a Pulse 
SEquence). This family is now seen to be just one mem- 
ber of a larger group of families. To choose between these 
it is necessary to look at higher order terms, and this is 
most conveniently achieved using the quaternion fidelity, 
Equation 0. As a baseline we take the fidelity of a sin- 
gle off-resonance 9<p pulse compared with its (ideal) on- 
resonance form, 



T « 1 + f 



(16) 



where terms in / 4 and higher have been neglected. Note 
that the fidelity only contains even order terms in / as 
the composite pulse perform symmetrically for positive 
and negative values of /. 

As expected all members of our general group of so- 
lutions result in much better fidelities; in particular the 
term in / 2 is always completely removed. The behaviour 
of the term in / 4 is much more complicated, but it can 
be shown that this term depends only on the value of 
n — n\ ~ n2 + ri3, that is the total number of additional 
2ir rotations performed by the composite pulse sequence, 
and has the smallest absolute value when the three in- 
tegers are chosen so that n = 0. As our previous values 
(m = 1, ri2 = 1 and ri^ = 0) are the smallest numbers 
that fit this criterion, it seems that the CORPSE family of 
pulse sequences is indeed the best member of this group. 
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TABLE I: Pulse rotation angles for a CORPSE composite pulse 
with a target rotation of 8 X ; CORPSE pulse phases are +x, —x, 
+x. 









FIG. 1: Fidelity of simple (dashed line) and CORPSE compos- 
ite pulses (solid line) as a function of the off-resonance fraction 
/ for pulses with a target rotation angle of 180°. 



The only other family of interest is that with n\ = 0, 
ri2 = 1 and U3 = 0, previously referred to as SHORT- 
CORPSE Q; while this performs less well than CORPSE 
it is somewhat shorter. Numerical values of individual 
pulse rotation angles for CORPSE sequences with a vari- 
ety of target angles are given in Table |. 

The performance of the CORPSE sequence for a f80° 
pulse is demonstrated in Figure Q the CORPSE pulse per- 
forms better than a simple pulse as long as |/| < 0.663. 
For smaller values of the effective range of / is reduced, 
but not dramatically so: for a 30° pulse the CORPSE pulse 
outperforms a simple pulse as long as |/| < 0.297. 



IV. PULSE LENGTH ERRORS 

A similar approach can be used to develop composite 
pulses to tackle pulse length errors. As before we begin 
with a sequence of three pulses, but the subsequent de- 
velopment is quite different. In particular we allow the 
three pulses to have arbitrary phase angles, as well as 
arbitrary rotation angles, thus giving us six variable pa- 
rameters, although this number is soon reduced to three. 

The quaternion corresponding to each pulse takes the 
simple form 



{cos(0'/2), sin(072){cos(0), sin(</>), 0}}. (17) 



where 0' = 0(1 + g) is the actual rotation angle achieved 
by a pulse with nominal rotation angle 0, and g is the 
fractional error in the pulse power. The quaternion de- 
scribing the composite pulse is very complicated, but can 
be simplified by restricting attention to the time symmet- 
ric case, where 0i = 03 and cj>i — 03. This automatically 
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ensures that the composite quaternion has no z compo- 
nent, as any time symmetric sequence of rotations about 
axes in the xy plane is itself a rotation about an axis in 
the xy plane. 

Even after this simplification, the composite quater- 
nion remains extremely complicated. To make further 
progress we note that a composite pulse of this kind has 
been previously described for the case of a 180° rotation: 
the sequence 



180 60 I8O300 180 



>G0 



(18) 



will perform a 180, T ro tation with compensation for pulse 
length errors (see p9j, but note the corrected phase an- 
gles). It seems likely that other members of this fam- 
ily will have either 9\ = tt or 62 = vr; both possibili- 
ties were initially explored, but the second choice seemed 
more productive and forms the basis of our subsequent 
work. 

As before, the composite quaternion can be expanded 
as a Maclaurin series in g, and it is most useful to con- 
centrate on the first order error term. This can be set 
equal to zero by choosing 



b\ ± arccos(— it/2Q\) 



(19) 



and, for consistency with equation |l8[ we will use the 
minus sign in future. Sequences obeying this equation 
will be insensitive to pulse length errors; the rotation and 
phase angle can then be adjusted by choosing suitable 
values for 9% and <p±. As before we will derive values for 
a 9 X pulse; pulses with other phase angles can be obtained 
by offsetting all the phase angles by the desired amount. 

Solving these equations is complex, but the solutions 
are fairly straightforward: 



h = "3 = arcsmc 



h = <p3 = arccos 



2cos(0/2) 



-TT COS Pi 



Wi sin(0/2) 
t>2 = <Pi — arccos(— 7r/20i) 



(20) 
(21) 
(22) 
(23) 



where sinc(x) is defined as sin(x)/x. We refer to this as 
a Short Composite Rotation For Undoing Length Over 
and Under Shoot or SCROFULOUS sequence. 

Numerical values of individual pulse rotation and 
phase angles for a variety of target angles are given in 
Table |fi|. The performance of SCROFULOUS and plain 
180° pulses are compared in Figure 0. 



V. THE BB1 FAMILY 

Another approach to composite pulse design has been 
described by Wimperis [p0| . While sequences such as 
corpse and scrofulous seek a single composite pulse 
which performs the desired rotation with reduced sensi- 
tivity to errors, an alternative approach is to combine a 
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TABLE II: Pulse rotation and phase angles for a SCROFULOUS 
composite pulse with a target rotation of 8 X ; note that 03 = 0i 
and 4>3 = (pi. 
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FIG. 2: Fidelity of simple (dashed line) and SCROFULOUS com- 
posite pulses (solid line) as a function of the fractional pulse 
length error g for pulses with a target rotation angle of 180° . 



naive pulse, which performs the desired rotation, with a 
sequence of error-correcting pulses, which partially com- 
pensate for imperfections. This approach appears to sim- 
plify the design of composite pulse sequences, and for the 
case of pulse length errors produces excellent results. 

The form of the error correcting pulse sequence is quite 
tightly constrained, as it must have no overall effect in the 
absence of errors, but it must also retain sufficient flexi- 
bility that it can act against errors when they do occur. 
Here we concentrate on one particular form suggested by 
Wimperis pi 



18O 01 36O 02 18O 01 



(24) 



where the values of <f>\ and <f>2 remain to be determined. 
When placed in front of a 9 X pulse Wimperis refers to 
the entire sequence as BB1 ^0|, but as we will generalise 
his approach we refer to the error correcting sequence 
(Equation ^J) as Wl. 

As before we evaluate the quaternion for the composite 
rotation (Wl followed by a X pulse) in the presence of 
pulse length errors, and then expand this quaternion as 
a Maclaurin series in g, the fractional error in the pulse 
power. The y and z components in the first order error 
term are easily removed by setting 4>2 — 30i ; the remain- 
ing components can then be eliminated by choosing 



± arccos 



47T 



(25) 



The positive solution is then identical to that previously 
described po| . Examining higher order error terms shows 
that this pulse sequence is even better than it first ap- 
pears, as these choices also completely remove the second 
order error terms. As discussed below, this effect appears 
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TABLE III: Pulse phase angles for a Wl correction sequence 
with a target rotation of 8 X ; pulse rotation angles are 9\ — 
180° and 62 = 360°. 
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FIG. 3: Fidelity of simple (dashed line) and BB1 composite 
pulses (solid line) as a function of the fractional pulse length 
error g for pulses with a target rotation angle of 180° . 



to be a property of the Wl sequence and its close rela- 
tions. 

It is easy to imagine a range of variations of the BB1 
sequence. Most simply the Wl error correction sequence 
can be placed after the 9 X pulse, instead of before it. 
Unsurprisingly this has no effect: the solution is the same 
as before, and the performance of this reversed sequence 
is identical to that of BB1. More surprisingly the Wl 
sequence can be place in the middle of the 9 X pulse, so 
that the overall sequence 



(0/2), Wl {0/2) x 



(26) 



is time symmetric. Indeed the Wl pulse can be placed 
at any point within the 6 X pulse, with almost identical 
effects. The form of the composite quaternion depends 
slightly on where the Wl pulse is placed, but the fidelity 
of the pulse sequence is unchanged: all error terms below 
sixth order are cancelled, with the size of the sixth order 
term depending on the value of 9. 

Numerical values of pulse phase angles for a variety of 
target angles are given in Table jTj. The performance of 
the BB1 and plain 180° pulses are compared in Figure ||. 
For all target angles below 180° the BB1 composite pulse 
outperforms a simple pulse when |<7| < 1. 

Another simple variation is to use two or more error 
correcting sequences; as before these can be placed at 
various different points around or within the 9 X pulse. 
For simplicity we assume that all the error correcting 
sequences are identical to one another, and have the same 
general form as Wl. In this case it can be shown that 
the correction sequences, which we call Wn, have phase 



angles given by 02 = 30i and 



± arccos 



4mt J 



(27) 



where n is the total number of sequences used. As before 
the fidelity is independent of where the Wn sequences are 
placed, but it does depend on the value of n. The second 
and fourth order error terms are cancelled in all cases, 
and the size of the sixth order error term now depends on 
both 9 and n. The smallest sixth order term is achieved 
when n = 2, but the term is not completely removed. 
The gain over n = 1 is fairly small, and in practice the 
simpler composite pulses based on the Wl sequence are 
likely to be the most effective. 

Having varied the position and number of the error 
correcting pulse sequences the next logical step is to vary 
their form. In principle any sequence that has no overall 
effect in the absence of errors could be used. In practice 
we find that many possible sequences allow the second 
order error term in the fidelity expression to be removed, 
but the simultaneous cancellation of second and fourth 
order errors seems to be a special feature of the Wn fam- 
ily of sequences. 

Given the success of this approach to tackling pulse 
length errors, it seems obvious to apply the method to 
tackle off-resonance effects. As yet, however, this ap- 
proach has had no success. 



VI. SIMULTANEOUS ERRORS 

So far we have only considered the case of either off- 
resonance effects or pulse length errors being present. In 
reality both problems may well occur simultaneously. It 
is therefore important to consider how such simultaneous 
errors might be tackled. Ideally we would like to design 
pulse sequences which can compensate for both problems 
at the same time; this, however, is a complicated and as 
yet unresolved problem, and here we simply analyse the 
sensitivity of each of our pulse sequences to the other 
kind of error. 

We proceed as before, calculating composite pulse and 
simple pulse quaternions in the presence of errors, and de- 
termining the quaternion fidelity. This fidelity can then 
be expanded as a Maclaurin series in the error, and the 
lower order terms examined. Note that this procedure 
still assumes that only one type of error is present at a 
time; in order to detail with the case where both errors 
are present simultaneously it would be possible to use a 
Maclaurin expansion in both errors, but this is unlikely 
to lead to much insight. Instead we will simply plot the 
fidelity as function of both errors for some chosen target 
angle. 

We begin by considering the response of the CORPSE 
pulse sequence to pulse length errors. In the absence of 
off-resonance effects the behaviour of CORPSE is trivial to 
calculate, as the three pulses are applied along the +x, 
—x and +x axes, so that the behaviour is identical to 
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that of a simple pulse. The behaviour of a 180° pulse in 
the presence of simultaneous errors is shown in Figure ^. 

The behaviour of the SCROFULOUS pulse sequence is 
difficult to calculate for general target rotation angles, 
due to the dependence of 6\ on the arcsinc function, and 
so we concentrate on the case of 180° pulses. For this 
case the dependence of the fidelity on off-resonance ef- 
fects is given by T « 1 — 2/ 2 , while a simple pulse has a 
fidelity JF « 1 - / 2 /2 (see Equation . In general SCRO- 
FULOUS is considerably more sensitive to off-resonance 
effects than plain pulses. 

Finally we consider the BB1 family of pulse sequences, 
taking the time-symmetrised version of BB1, Equa- 
tion |26|, as our standard. In this case we can solve the 
problem for any target rotation angle, and up to second 
order the result is identical to that of a plain pulse, Equa- 
tion [16]. Thus, unlike SCROFULOUS, the BB1 sequence 
achieves its impressive tolerance to pulse length errors 
at little or no cost in sensitivity to off-resonance effects. 
This is confirmed for simultaneous errors by Figure 0. 



implementation of a quantum computer must be con- 
cerned on some level with rotations on the Bloch sphere, 
and so composite pulse techniques may find very broad 
application in quantum computing. Composite pulses 
are not, however, a panacea, and some caution must be 
exercised in their use. 

The CORPSE pulse sequence appears to be the best 
approach for tackling small off-resonance errors (for large 
known off-resonance effects the resonance offset tailored, 
or rotten, scheme is preferable). For pulse length 
errors variations on the BB1 scheme of Wimperis |2C(| 
give the best results; the SCROFULOUS family of pulses 
is less effective, but does have the advantage of being 
considerably shorter. 
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